Exploring the diversity and evolutionary strategies of prophages in Hyphomicrobiales, comparing animal-associated with non-animal-associated bacteria.

The Hyphomicrobiales bacterial order (previously Rhizobiales) exhibits a wide range of lifestyle characteristics, including free-living, plant-association, nitrogen-fixing, and association with animals (Bartonella and Brucella). This study explores the diversity and evolutionary strategies of bacteriophages within the Hyphomicrobiales order, comparing animal-associated (AAB) with non-animal-associated bacteria (NAAB). We curated 560 high-quality complete genomes of 58 genera from this order and used the PHASTER server for prophage annotation and classification. For 19 genera with representative genomes, we curated 96 genomes and used the Defense-Finder server to summarize the type of anti-phage systems (APS) found in this order. We analyzed the genetic repertoire and length distributions of prophages, estimating evolutionary rates and comparing intact, questionable, and incomplete prophages in both groups. Analyses of best-fit parameters and bootstrap sensitivity were used to understand the evolutionary processes driving prophage gene content. A total of 1860 prophages distributed in Hyphomicrobiales were found, 695 in AAB and 1165 in the NAAB genera. The results revealed a similar number of prophages per genome in AAB and NAAB and a similar length distribution, suggesting shared mechanisms of genetic acquisition of prophage genes. Changes in the frequency of specific gene classes were observed between incomplete and intact prophages, indicating preferential loss or enrichment in both groups. The analysis of best-fit parameters and bootstrap sensitivity tests indicated a higher selection coefficient, induction rate, and turnover in NAAB genomes. We found 68 types of APS in Hyphomicrobiales; restriction modification (RM) and abortive infection (Abi) were the most frequent APS found for all Hyphomicrobiales, and within the AAB group. This classification of APS showed that NAAB genomes have a greater diversity of defense systems compared to AAB, which could be related to the higher rates of prophage induction and turnover in the latter group. Our study provides insights into the distributions of both prophages and APS in Hyphomicrobiales genomes, demonstrating that NAAB carry more defense systems against phages, while AAB show increased prophage stability and an increased number of incomplete prophages. These results suggest a greater role for domesticated prophages within animal-associated bacteria in Hyphomicrobiales.

Integrative and Conjugative Elements and Prophage DNA as Carriers of Resistance Genes in Erysipelothrix rhusiopathiae Strains from Domestic Geese in Poland.

Goose erysipelas is a serious problem in waterfowl breeding in Poland. However, knowledge of the characteristics of Erysipelothrix rhusiopathiae strains causing this disease is limited. In this study, the antimicrobial susceptibility and serotypes of four E. rhusiopathiae strains from domestic geese were determined, and their whole-genome sequences (WGSs) were analyzed to detect resistance genes, integrative and conjugative elements (ICEs), and prophage DNA. Sequence type and the presence of resistance genes and transposons were compared with 363 publicly available E. rhusiopathiae strains, as well as 13 strains of other Erysipelothrix species. Four strains tested represented serotypes 2 and 5 and the MLST groups ST 4, 32, 242, and 243. Their assembled circular genomes ranged from 1.8 to 1.9 kb with a GC content of 36-37%; a small plasmid was detected in strain 1023. Strains 1023 and 267 were multidrug-resistant. The resistance genes detected in the genome of strain 1023 were erm47, tetM, and lsaE-lnuB-ant(6)-Ia-spw cluster, while strain 267 contained the tetM and ermB genes. Mutations in the gyrA gene were detected in both strains. The tetM gene was embedded in a Tn916-like transposon, which in strain 1023, together with the other resistance genes, was located on a large integrative and conjugative-like element of 130 kb designated as ICEEr1023. A minor integrative element of 74 kb was identified in strain 1012 (ICEEr1012). This work contributes to knowledge about the characteristics of E. rhusiopathiae bacteria and, for the first time, reveals the occurrence of erm47 and ermB resistance genes in strains of this species. Phage infection appears to be responsible for the introduction of the ermB gene into the genome of strain 267, while ICEs most likely play a key role in the spread of the other resistance genes identified in E. rhusiopathiae.

Slings and arrows: sRNAs mediate intragenomic competition.

Bacterial genomes are littered with exogenous: competing DNA elements. Here, Sprenger et al. demonstrate that the Vibrio cholerae prophage VP882 modulates host functions via production of regulatory sRNAs to promote phage development. Alternatively, host sRNAs inhibit the VP882 lytic phase by specifically regulating phage genes.

Genetic recombination-mediated evolutionary interactions between phages of potential industrial importance and prophages of their hosts within or across the domains of Escherichia, Listeria, Salmonella, Campylobacter, and Staphylococcus.

BACKGROUND: The in-depth understanding of the role of lateral genetic transfer (LGT) in phage-prophage interactions is essential to rationalizing phage applications for human and animal therapy, as well as for food and environmental safety. This in silico study aimed to detect LGT between phages of potential industrial importance and their hosts. METHODS: A large array of genetic recombination detection algorithms, implemented in SplitsTree and RDP4, was applied to detect LGT between various Escherichia, Listeria, Salmonella, Campylobacter, Staphylococcus, Pseudomonas, and Vibrio phages and their hosts. PHASTER and RAST were employed respectively to identify prophages across the host genome and to annotate LGT-affected genes with unknown functions. PhageAI was used to gain deeper insights into the life cycle history of recombined phages. RESULTS: The split decomposition inferences (bootstrap values: 91.3-100; fit: 91.433-100), coupled with the Phi (0.0-2.836E-12) and RDP4 (P being well below 0.05) statistics, provided strong evidence for LGT between certain Escherichia, Listeria, Salmonella, and Campylobacter virulent phages and prophages of their hosts. The LGT events entailed mainly the phage genes encoding for hypothetical proteins, while some of these genetic loci appeared to have been affected even by intergeneric recombination in specific E. coli and S. enterica virulent phages when interacting with their host prophages. Moreover, it is shown that certain L. monocytogenes virulent phages could serve at least as the donors of the gene loci, involved in encoding for the basal promoter specificity factor, for L. monocytogenes. In contrast, the large genetic clusters were determined to have been simultaneously exchanged by many S. aureus prophages and some Staphylococcus temperate phages proposed earlier as potential therapeutic candidates (in their native or modified state). The above genetic clusters were found to encompass multiple genes encoding for various proteins, such as e.g., phage tail proteins, the capsid and scaffold proteins, holins, and transcriptional terminator proteins. CONCLUSIONS: It is suggested that phage-prophage interactions, mediated by LGT (including intergeneric recombination), can have a far-reaching impact on the co-evolutionary trajectories of industrial phages and their hosts especially when excessively present across microbially rich environments.

The complete genome sequence of unculturable Mycoplasma faucium obtained through clinical metagenomic next-generation sequencing.

Introduction: Diagnosing Mycoplasma faucium poses challenges, and it's unclear if its rare isolation is due to infrequent occurrence or its fastidious nutritional requirements. Methods: This study analyzes the complete genome sequence of M. faucium, obtained directly from the pus of a sternum infection in a lung transplant patient using metagenomic sequencing. Results: Genome analysis revealed limited therapeutic options for the M. faucium infection, primarily susceptibility to tetracyclines. Three classes of mobile genetic elements were identified: two new insertion sequences, a new prophage (phiUMCG-1), and a species-specific variant of a mycoplasma integrative and conjugative element (MICE). Additionally, a Type I Restriction-Modification system was identified, featuring 5'-terminally truncated hsdS pseudogenes with overlapping repeats, indicating the potential for forming alternative hsdS variants through recombination. Conclusion: This study represents the first-ever acquisition of a complete circularized bacterial genome directly from a patient sample obtained from invasive infection of a primary sterile site using culture-independent, PCR-free clinical metagenomics.

Bacteriophages avoid autoimmunity from cognate immune systems as an intrinsic part of their life cycles.

Dormant prophages protect lysogenic cells by expressing diverse immune systems, which must avoid targeting their cognate prophages upon activation. Here we report that multiple Staphylococcus aureus prophages encode Tha (tail-activated, HEPN (higher eukaryotes and prokaryotes nucleotide-binding) domain-containing anti-phage system), a defence system activated by structural tail proteins of incoming phages. We demonstrate the function of two Tha systems, Tha-1 and Tha-2, activated by distinct tail proteins. Interestingly, Tha systems can also block reproduction of the induced tha-positive prophages. To prevent autoimmunity after prophage induction, these systems are inhibited by the product of a small overlapping antisense gene previously believed to encode an excisionase. This genetic organization, conserved in S. aureus prophages, allows Tha systems to protect prophages and their bacterial hosts against phage predation and to be turned off during prophage induction, balancing immunity and autoimmunity. Our results show that the fine regulation of these processes is essential for the correct development of prophages' life cycle.

Phage biology: The ins and outs of prophages in bacterial populations.

Bacterial genomes often harbor integrated viruses (prophages), which provide novel functions but also lyse cells under stressful conditions. A new paper combines mathematical models with experimental evolution to determine how prophages are maintained in bacterial populations despite their fitness costs.

Prophage maintenance is determined by environment-dependent selective sweeps rather than mutational availability.

Prophages, viral sequences integrated into bacterial genomes, can be beneficial and costly. Despite the risk of prophage activation and subsequent bacterial death, active prophages are present in most bacterial genomes. However, our understanding of the selective forces that maintain prophages in bacterial populations is limited. Combining experimental evolution with stochastic modeling, we show that prophage maintenance and loss are primarily determined by environmental conditions that alter the net fitness effect of a prophage on its bacterial host. When prophages are too costly, they are rapidly lost through environment-specific sequences of selective sweeps. Conflicting selection pressures that select against the prophage but for a prophage-encoded accessory gene can maintain prophages. The dynamics of prophage maintenance additionally depend on the sociality of this accessory gene. Prophage-encoded genes that exclusively benefit the lysogen maintain prophages at higher frequencies compared with genes that benefit the entire population. That is because the latter can protect phage-free "cheaters," reducing the benefit of maintaining the prophage. Our simulations suggest that environmental variation plays a larger role than mutation rates in determining prophage maintenance. These findings highlight the complexity of selection pressures that act on mobile genetic elements and challenge our understanding of the role of environmental factors relative to random chance events in shaping the evolutionary trajectory of bacterial populations. By shedding light on the key factors that shape microbial populations in the face of environmental changes, our study significantly advances our understanding of the complex dynamics of microbial evolution and diversification.

Genome Analysis of Epsilon CrAss-like Phages.

CrAss-like phages play an important role in maintaining ecological balance in the human intestinal microbiome. However, their genetic diversity and lifestyle are still insufficiently studied. In this study, a novel CrAssE-Sib phage genome belonging to the epsilon crAss-like phage genomes was found. Comparative analysis indicated that epsilon crAss-like phages are divided into two putative genera, which were proposed to be named Epsilonunovirus and Epsilonduovirus; CrAssE-Sib belongs to the former. The crAssE-Sib genome contains a diversity-generating retroelement (DGR) cassette with all essential elements, including the reverse transcriptase (RT) and receptor binding protein (RBP) genes. However, this RT contains the GxxxSP motif in its fourth domain instead of the usual GxxxSQ motif found in all known phage and bacterial DGRs. RBP encoded by CrAssE-Sib and other Epsilonunoviruses has an unusual structure, and no similar phage proteins were found. In addition, crAssE-Sib and other Epsilonunoviruses encode conserved prophage repressor and anti-repressors that could be involved in lysogenic-to-lytic cycle switches. Notably, DNA primase sequences of epsilon crAss-like phages are not included in the monophyletic group formed by the DNA primases of all other crAss-like phages. Therefore, epsilon crAss-like phage substantially differ from other crAss-like phages, indicating the need to classify these phages into a separate family.

Assessing phage-host population dynamics by reintroducing virulent viruses to synthetic microbiomes.

Microbiomes feature complex interactions between diverse bacteria and bacteriophages. Synthetic microbiomes offer a powerful way to study these interactions; however, a major challenge is obtaining a representative bacteriophage population during the bacterial isolation process. We demonstrate that colony isolation reliably excludes virulent viruses from sample sources with low virion-to-bacteria ratios such as feces, creating "virulent virus-free" controls. When the virulent dsDNA virome is reintroduced to a 73-strain synthetic gut microbiome in a bioreactor model of the human colon, virulent viruses target susceptible strains without significantly altering community structure or metabolism. In addition, we detected signals of prophage induction that associate with virulent predation. Overall, our findings indicate that dilution-based isolation methods generate synthetic gut microbiomes that are heavily depleted, if not devoid, of virulent viruses and that such viruses, if reintroduced, have a targeted effect on community assembly, metabolism, and prophage replication.

Prophage terminase with tRNase activity sensitizes Salmonella enterica to oxidative stress.

Phage viruses shape the evolution and virulence of their bacterial hosts. The Salmonella enterica genome encodes several stress-inducible prophages. The Gifsy-1 prophage terminase protein, whose canonical function is to process phage DNA for packaging in the virus head, unexpectedly acts as a transfer ribonuclease (tRNase) under oxidative stress, cleaving the anticodon loop of tRNALeu. The ensuing RNA fragmentation compromises bacterial translation, intracellular survival, and recovery from oxidative stress in the vertebrate host. S. enterica adapts to this transfer RNA (tRNA) fragmentation by transcribing the RNA repair Rtc system. The counterintuitive translational arrest provided by tRNA cleavage may subvert prophage mobilization and give the host an opportunity for repair as a way of maintaining bacterial genome integrity and ultimately survival in animals.

CRISPRpi: Inducing and Curing Prophage Using the CRISPR Interference.

We present here a CRISPR-interference-based protocol to trigger prophage induction, even for non-inducible prophages. This method can also be used to cure the prophage from the bacterial host. The method is based on silencing of the phage's repressor transcription, thanks to CRISPR interference. Plasmid electroporation is used to bring the CRISPRi system into the bacteria, specifically on a plasmid carrying spacers targeting the prophage repressor. This method enables prophage induction and curation in a week or two with a high efficiency.

The spontaneously produced lysogenic prophage phi456 promotes bacterial resistance to adverse environments and enhances the colonization ability of avian pathogenic Escherichia coli strain DE456.

In the last decade, prophages that possess the ability of lysogenic transformation have become increasingly significant. Their transfer and subsequent activity in the host have a significant impact on the evolution of bacteria. Here, we investigate the role of prophage phi456 with high spontaneous induction in the bacterial genome of Avian pathogenic Escherichia coli (APEC) DE456. The phage particles, phi456, that were released from DE456 were isolated, purified, and sequenced. Additionally, phage particles were no longer observed either during normal growth or induced by nalidixic acid in DE456Δphi456. This indicated that the released phage particles from DE456 were only phi456. We demonstrated that phi456 contributed to biofilm formation through spontaneous induction of the accompanying increase in the eDNA content. The survival ability of DE456Δphi456 was decreased in avian macrophage HD11 under oxidative stress and acidic conditions. This is likely due to a decrease in the transcription levels of three crucial genes-rpoS, katE, and oxyR-which are needed to help the bacteria adapt to and survive in adverse environments. It has been observed through animal experiments that the presence of phi456 in the DE456 genome enhances colonization ability in vivo. Additionally, the number of type I fimbriae in DE456Δphi456 was observed to be reduced under transmission electron microscopy when compared to the wild-type strain. The qRT-PCR results indicated that the expression levels of the subunit of I fimbriae (fimA) and its apical adhesin (fimH) were significantly lower in DE456Δphi456. Therefore, it can be concluded that phi456 plays a crucial role in helping bacterial hosts survive in unfavorable conditions and enhancing the colonization ability in DE456.

Diverse Prophage Elements of Salmonella enterica Serovars Show Potential Roles in Bacterial Pathogenicity.

Nontyphoidal salmonellosis is an important foodborne and zoonotic infection that causes significant global public health concern. Diverse serovars are multidrug-resistant and encode several virulence indicators; however, little is known on the role prophages play in driving these traits. Here, we extracted prophages from seventy-five Salmonella genomes which represent the fifteen important serovars in the United Kingdom. We analyzed the intact prophages for the presence of virulence genes and established their genomic relationships. We identified 615 prophages from the Salmonella strains, from which 195 prophages are intact, 332 are incomplete, while 88 are questionable. The average prophage carriage was found to be 'extreme' in S. Heidelberg, S. Inverness, and S. Newport (10.2-11.6 prophages/strain), 'high' in S. Infantis, S. Stanley, S. Typhimurium, and S. Virchow (8.2-9.0 prophages/strain), 'moderate' in S. Agona, S. Braenderup, S. Bovismorbificans, S. Choleraesuis, S. Dublin, and S. Java (6.0-7.8 prophages/strain), and 'low' in S. Javiana and S. Enteritidis (5.8 prophages/strain). Cumulatively, 61 virulence genes (1500 gene copies) were detected from representative intact prophages and linked to Salmonella delivery/secretion system (42.62%), adherence (32.7%), magnesium uptake (3.88%), regulation (5%), stress/survival (1.6%), toxins (10%), and antivirulence (1.6%). Diverse clusters were formed among the intact prophages and with bacteriophages of other enterobacteria, suggesting different lineages and associations. Our work provides a strong body of data to support the contributions diverse prophages make to the pathogenicity of Salmonella, including thirteen previously unexplored serovars.

Harnessing the Diversity of Burkholderia spp. Prophages for Therapeutic Potential.

Burkholderia spp. are often resistant to antibiotics, and infections with these organisms are difficult to treat. A potential alternative treatment for Burkholderia spp. infections is bacteriophage (phage) therapy; however, it can be difficult to locate phages that target these bacteria. Prophages incorporated into the bacterial genome have been identified within Burkholderia spp. and may represent a source of useful phages for therapy. Here, we investigate whether prophages within Burkholderia spp. clinical isolates can kill conspecific and heterospecific isolates. Thirty-two Burkholderia spp. isolates were induced for prophage release, and harvested phages were tested for lytic activity against the same 32 isolates. Temperate phages were passaged and their host ranges were determined, resulting in four unique phages of prophage origin that showed different ranges of lytic activity. We also analyzed the prophage content of 35 Burkholderia spp. clinical isolate genomes and identified several prophages present in the genomes of multiple isolates of the same species. Finally, we observed that Burkholdera cenocepacia isolates were more phage-susceptible than Burkholderia multivorans isolates. Overall, our findings suggest that prophages present within Burkholderia spp. genomes are a potentially useful starting point for the isolation and development of novel phages for use in phage therapy.

Re-analysis of an outbreak of Shiga toxin-producing Escherichia coli O157:H7 associated with raw drinking milk using Nanopore sequencing.

The aim of this study was to compare Illumina and Oxford Nanopore Technology (ONT) sequencing data to quantify genetic variation to assess within-outbreak strain relatedness and characterise microevolutionary events in the accessory genomes of a cluster of 23 genetically and epidemiologically linked isolates related to an outbreak of Shiga toxin-producing Escherichia coli O157:H7 caused by the consumption of raw drinking milk. There were seven discrepant variants called between the two technologies, five were false-negative or false-positive variants in the Illumina data and two were false-negative calls in ONT data. After masking horizontally acquired sequences such as prophages, analysis of both short and long-read sequences revealed the 20 isolates linked to the outbreak in 2017 had a maximum SNP distance of one SNP between each other, and a maximum of five SNPs when including three additional strains identified in 2019. Analysis of the ONT data revealed a 47 kbp deletion event in a terminal compound prophage within one sample relative to the remaining samples, and a 0.65 Mbp large chromosomal rearrangement (inversion), within one sample relative to the remaining samples. Furthermore, we detected two bacteriophages encoding the highly pathogenic Shiga toxin (Stx) subtype, Stx2a. One was typical of Stx2a-phage in this sub-lineage (Ic), the other was atypical and inserted into a site usually occupied by Stx2c-encoding phage. Finally, we observed an increase in the size of the pO157 IncFIB plasmid (1.6 kbp) in isolates from 2019 compared to those from 2017, due to the duplication of insertion elements within the plasmids from the more recently isolated strains. The ability to characterize the accessory genome in this way is the first step to understanding the significance of these microevolutionary events and their impact on the genome plasticity and virulence between strains of this zoonotic, foodborne pathogen.

Prophage proteins alter long noncoding RNA and DNA of developing sperm to induce a paternal-effect lethality.

The extent to which prophage proteins interact with eukaryotic macromolecules is largely unknown. In this work, we show that cytoplasmic incompatibility factor A (CifA) and B (CifB) proteins, encoded by prophage WO of the endosymbiont Wolbachia, alter long noncoding RNA (lncRNA) and DNA during Drosophila sperm development to establish a paternal-effect embryonic lethality known as cytoplasmic incompatibility (CI). CifA is a ribonuclease (RNase) that depletes a spermatocyte lncRNA important for the histone-to-protamine transition of spermiogenesis. Both CifA and CifB are deoxyribonucleases (DNases) that elevate DNA damage in late spermiogenesis. lncRNA knockdown enhances CI, and mutagenesis links lncRNA depletion and subsequent sperm chromatin integrity changes to embryonic DNA damage and CI. Hence, prophage proteins interact with eukaryotic macromolecules during gametogenesis to create a symbiosis that is fundamental to insect evolution and vector control.

Speeding genomic island discovery through systematic design of reference database composition.

BACKGROUND: Genomic islands (GIs) are mobile genetic elements that integrate site-specifically into bacterial chromosomes, bearing genes that affect phenotypes such as pathogenicity and metabolism. GIs typically occur sporadically among related bacterial strains, enabling comparative genomic approaches to GI identification. For a candidate GI in a query genome, the number of reference genomes with a precise deletion of the GI serves as a support value for the GI. Our comparative software for GI identification was slowed by our original use of large reference genome databases (DBs). Here we explore smaller species-focused DBs. RESULTS: With increasing DB size, recovery of our reliable prophage GI calls reached a plateau, while recovery of less reliable GI calls (FPs) increased rapidly as DB sizes exceeded ~500 genomes; i.e., overlarge DBs can increase FP rates. Paradoxically, relative to prophages, FPs were both more frequently supported only by genomes outside the species and more frequently supported only by genomes inside the species; this may be due to their generally lower support values. Setting a DB size limit for our SMAll Ranked Tailored (SMART) DB design speeded runtime ~65-fold. Strictly intra-species DBs would tend to lower yields of prophages for small species (with few genomes available); simulations with large species showed that this could be partially overcome by reaching outside the species to closely related taxa, without an FP burden. Employing such taxonomic outreach in DB design generated redundancy in the DB set; as few as 2984 DBs were needed to cover all 47894 prokaryotic species. CONCLUSIONS: Runtime decreased dramatically with SMART DB design, with only minor losses of prophages. We also describe potential utility in other comparative genomics projects.

Escherichia coli CRISPR arrays from early life fecal samples preferentially target prophages.

CRISPR-Cas systems are defense mechanisms against phages and other nucleic acids that invade bacteria and archaea. In Escherichia coli, it is generally accepted that CRISPR-Cas systems are inactive in laboratory conditions due to a transcriptional repressor. In natural isolates, it has been shown that CRISPR arrays remain stable over the years and that most spacer targets (protospacers) remain unknown. Here, we re-examine CRISPR arrays in natural E. coli isolates and investigate viral and bacterial genomes for spacer targets using a bioinformatics approach coupled to a unique biological dataset. We first sequenced the CRISPR1 array of 1769 E. coli isolates from the fecal samples of 639 children obtained during their first year of life. We built a network with edges between isolates that reflect the number of shared spacers. The isolates grouped into 34 modules. A search for matching spacers in bacterial genomes showed that E. coli spacers almost exclusively target prophages. While we found instances of self-targeting spacers, those involving a prophage and a spacer within the same bacterial genome were rare. The extensive search for matching spacers also expanded the library of known E. coli protospacers to 60%. Altogether, these results favor the concept that E. coli's CRISPR-Cas is an antiprophage system and highlight the importance of reconsidering the criteria use to deem CRISPR-Cas systems active.

Oral bacteriophages: metagenomic clues to interpret microbiomes.

Bacteriophages are bacterial viruses that are distributed throughout the environment. Lytic phages and prophages in saliva, oral mucosa, and dental plaque interact with the oral microbiota and can change biofilm formation. The interactions between phages and bacteria can be considered a portion of oral metagenomics. The metagenomic profile of the oral microbiome indicates various bacteria. Indeed, there are various phages against these bacteria in the oral cavity. However, some other phages, like phages against Absconditabacteria, Chlamydiae, or Chloroflexi, have not been identified in the oral cavity. This review gives an overview of oral bacteriophage and used for metagenomics. Metagenomics of these phages deals with multi-drug-resistant bacterial plaques (biofilms) in oral cavities and oral infection. Hence, dentists and pharmacologists should know this metagenomic profile to cope with predental and dental infectious diseases.